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We investigate the stress-driven morphological instability of epitaxially growing multilayer films, 
which are coherent and dislocation- free. We construct a direct elastic analysis, from which we 
^~>) ' determine the elastic state of the system recursively in terms of that of the old states of the buried 

' layers. In turn, we use the result for the elastic state to derive the morphological evolution equation 

, of surface profile to first order of perturbations, with the solution explicitly expressed by the growth 

K ■ conditions and material parameters of all the deposited layers. We apply these results to two kinds 

of multilayer structures. One is the alternating tensile/compressive multilayer structure, for which 
we determine the effective stability properties, including the effect of varying surface mobility in 
different layers, its interplay with the global misfit of the multilayer film, and the influence of 
asymmetric structure of compressive and tensile layers on the system stability. The nature of the 
' (— I ' I asymmetry properties found in stability diagrams is in agreement with experimental observations. 

. The other multilayer structure that we study is one composed of stacked strained/spacer layers. We 

' also calculate the kinetic critical thickness for the onset of morphological instability and obtain its 

C ' reduction and saturation as number of deposited layers increases, which is consistent with recent 

experimental results. Compared to the single-layer film growth, the behavior of kinetic critical 
^ thickness shows deviations for upper strained layers. 
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I. INTRODUCTION 



In recent years, much attention has been focused on the heteroepitaxial growth of multilayer films for the purpose of 
'■^ fabricating high-quality low-dimensional structures.^ Currently, two kinds of multilayers are being actively studied: (i) 
^ ' the alternating tensile/compressive multilayers exploited as multiple quantum wells or short-period superlattices;^^* 
O I and (ii) multisheets of coherently strained layers (or islands) alternately separated by spacer layers of substrate 
^ ■ materials (strained/spacer arrays). ^'^"^^ Experimentally, these multilayer structures are grown coherently on the sub- 
strates, without the generation of undesirable misfit dislocations. When instability occurs and develops during the 

■ multilayer films growth, one of the most important phenomena is the formation of self-organized nanostructures. 
^ ' For the tensile/compressive multilayers, lateral modulations have been found during the epitaxial growth, show- 
ing as the modulated pattern of alternating vertical or columnar structures with strong contrast compositions. ■^^^ 
These modulation structures that form spontaneously are very useful for the quantum wires fabrication. For the 

[ strained/spacer multilayer structure, a remarkable phenomenon is the vertical self-organization of multisheet arrays 

■ of three-dimensional (3D) islands (quantum dots).^'^"^^ The islands are vertically correlated in successive strained 
CN| [ layers, with the increasing uniformity (including island size and density) as the number of layers increases. These 

. low-dimensional nanostructures (such as quantum wires and quantum dots), realized through the development of 
multilayer film instability, exhibit novel electronic and optical properties, including the band structure, carrier trans- 
port, etc., and are of intense interests both for fundamental studies and for application to devices with enhanced 
optoelectronic performance. 

The mechanism of occurrence of instability in strained-layer growth is believed to be stress-driven. When a single- 
^ ' layer strained film (with single component) is grown epitaxially on a substrate, the morphological instability develops 
Q , due to the nonuniform strain energy distribution with respect to any surface corrugations and the surface diffusion 
O ■ processes driven by the difference of local strain energy.^^'^° When other layers composed of different materials, and 
^ ' consequently with different misfit strains with respect to the substrate, are deposited subsequently, as in the multilayer 
growth, the coupling of strain fields caused by buried layers and the growing one determines the morphological 
evolution of surface profile. This coupling is the source for self-organized growth of vertical islands array in the 
^ , strained/spacer multilayer. ^'^^ For the tensile/compressive multilayer films, the phenomenon of lateral composition 
modulations can be explained by the morphological instability caused by the coupling strains, leading to the surface 
undulations and then to both the layer thickness modulations and composition modulations in the overall film through 
the stacking of layers. ^'^ 

Although in the coherent multilayer systems the processes of morphological evolution are mainly dominated by 
elastic effects, our understanding of the mechanisms is incomplete and not well developed due to the complexity of 
the multilayer structures. One of the main difficulties of the theoretical investigations is in determining the elastic 
fields of the system, especially the surface strain field caused by all the buried layers or islands. First, the system 
is inhomogeneous due to the stacking of layers with different materials; second, the geometry of this system, which 
is composed of top surface, finite thickness of film layers and the semi-infinite substrate, is not symmetric and is 
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therefore different from the case of bulk materials. Moreover, the interfaces between different layers arc generally not 
planar (except for the interfaces between strained and thick enough spacer layers), and then the system cannot be 
simplified to an equivalent and solvable semi-infinite solid as in the study of single-layer films. Thus, it is very 
difficult to exactly calculate the elastic fields of the system. 

The theoretical work for multilayer systems is rather sparse, and some approximations have to be used. For the 
tensile/compressive multilayer structure, Sridhar et al.^^ considered an infinite periodic and symmetric system with 
the dynamics controlled by interfacial diffusion, and focused on the interior layers that are far away from both surface 
and substrate, instead of the surface morphology that we are interested in this paper. In the recent work of Shilkrot et 
al.,'^'^ the surface diffusion and morphological evolution is discussed, with the determination of stability for multilayer 
films. The Eshelby procedure^^ is used for each buried interface, which is assumed to separately contribute to surface 
strain. Then the problem is approximated to be that of a misfit inclusion in a semi-infinite homogeneous isotropic 
elastic medium with the presence of a free surface. For a system of stacked coherent strained islands (strained/spacer 
multilayer), in the work of Tersoff et al.^ each buried island is approximately treated as a point defect, and its 
contribution to surface strain is determined by using the continuum elasticity solutions of Maradudin and Wallis.^^ 
Their results explain the vertical correlation between islands of successive layers, showing the self-organization of the 
system. The size and shape of stacked islands also exhibit the self-organized effects, as studied in a later work by Liu 
et al.}'^ 

Here, we propose a recursive procedure and develop a continuum, dynamical model to directly describe the mor- 
phological evolution of these coherent multilayer structures during the epitaxial growth. Using the isotropic linear 
elasticity theory, we calculate in section II the clastic fields of the multilayer system directly from the linearized 
boundary conditions at all the buried interfaces and the top surface. Based on our elastic analysis, the coupling and 
influence of all the buried layers and interfaces are determined in section III to describe the evolution of surface mor- 
phological perturbation. Within this framework, we investigate the stability properties of multilayer growing films by 
performing a linear analysis. Also, the effects of different material parameters and growth conditions on morphological 
stability/instability of multilayer films are determined and discussed. 

We first apply in section IV these results to the alternating tensile/compressive multilayer system, determining 
the effective stability properties. We investigate the effect of different surface mobilities of strained layers on the 
stability results, and the interplay between the relative mobility and the global misfit of multilayer film (especially 
the asymmetry between globally tensile and compressive films). In the absence of global strain, i.e., with a global 
strain-balanced condition, we find an interesting property that the system is more stable for the asymmetric structure 
(with respect to layer thickness and misfit) of compressive and tensile layers. The effects of relative surface tension 
and growth temperature are also studied. 

Then in section V we investigate the strained/spacer multilayer structure (for the case of flat surface of spacer layers), 
and calculate the kinetic critical thickness for the onset of 3D instability. We obtain a reduction and saturation of 
critical thickness in the upper strained layers, well reproducing the recent experimental observations. We also find 
that the behavior of kinetic critical thickness for upper strained layers deviates from that of single-layer film growth 
for large misfit and small deposition rate. 

II. MODEL AND ELASTIC FORMULATION 

Assume that a multilayer film is coherently deposited on a semi-infinite substrate (denoted by index 0) and consists 
of different layers i = 1, 2, k. In the general case, material of the deposited layer i can be different from that of 
the substrate and that of all the other layers, and its misfit with respect to the substrate is denoted as Cj. As we know 
from the hctcroepitaxial growth of a single- layer hlm,^''^^^'^'''^^ the nonzero misfit strain can be relieved through the 
development of morphological instability, leading to a nonplanar surface profile z = hi(x, y) = ^i+X]q ^il*?) G'^\>{iqxX+ 
iqyV) for the first layer 1, with a coherent interface at z = ho = C,{x,y) = C(q) 6xp(iga;a; -|- iqyy) (the average 

position C = 0) between film and substrate, where l\ is the average thickness of first layer and hi {() denotes the 
surface (interface) morphological perturbation. Here we only consider the morphological modulations of the film, 
corresponding to the deposition of only a single component in each layer, and do not consider the compositional 
evolution which is important in alloy growth systems and which results in more complicated properties due to its 
coupling with the morphological evolution. 

When the second and higher layers are deposited with different materials and thus with different misfits, the 
undulating surface is buried as a frozen rough interface if the growth temperature T is not too high making interlayer 
diffusion negligible. Then the elastic state of the system is altered as a result of new deposition, and the morphological 
evolution of the new free surface is different from that of the first single layer due to the coupling of strain fields in 
different layers and the influence of buried interfaces on the elastic field and energy of the multilayer film as a whole. For 
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a fc-layer system, the growing surface profile is represented as hk{x, y, t) = Y^^=i h + Vkt + J^q ^fe(q, t) G'y^v(ilxX + iqyy), 

with li {\ < i < k) the average thickness of ith layer underneath, Vk the deposition rate of top layer k, and hk the 
top surface perturbation. When we consider the atomic diffusion along the top surface and the uniform deposition 
of material as the dominant mass-transport processes, the evolution equation of the surface morphological profile is 
given by 

where = 1 + |V/ifep denotes the determinant of surface metric for top layer fc, and = Dg^kNg^k/ksTN^ j. is 
the surface mobility with Dg^k the surface diffusivity, ks the Boltzmann constant. Ns^k and Ny^k respectively the 
number density of atoms per unit surface area and per unit volume of layer k. Both the conserved thermal noise t]c 
and nonconserved deposition noise r]k are included in the above evolution equation, but are neglected in the following 
studies since they do not alter the results of linear analysis presented here. In Eq. (1), the free energy of the system 
consists of two terms: 



= 7fc / d\^+ [ " d\£. (2) 

J J — oo 



The first term is the surface energy with jk the surface tension at the top free surface of layer k, and tends to stabilize 
the film morphology. £ in the second term represents the strain energy density, and is, in linear elasticity theory, 
£ = Safi^pcra^f^ipf^ {ct, P = 2;, V, z), where aa/s is the stress tensor. For an elastically isotropic system that 

we assume here, the elastic compliance tensor Sap^p = Sa^Si3p{l + i')/E — Sa/sS^pv/E with Young's modulus E and 
Poisson ratio v. The form of total free energy in Eq. (2) implies that wc have neglected the wetting effect of each 
deposited layer, which is not dominant when the layer thickness is assumed to be large enough. The wetting effect is 
related to the interaction between two successive layers or between film and substrate, and is important for thin films 
with thickness less than 3 ML.^® 

Substituting Eq. (2) into Eq. (1), wc obtain the expression for surface profile evolution of top growing layer k: 

dhk/dt ^Vk^/gk^l[lkH + £k{hk,hk-i, ■■■X)] +Vk, (3) 

where k is the surface curvature, and £k is the elastic energy density evaluated at surface z = hk and is dependent on 
all the surface and interface morphologies due to the coupling of strain fields. Since here we focus on the stability of 
the multilayer system and hence the early evolution regime with small perturbations, to determine the corresponding 
important properties we only need the first-order evolution equation for surface morphology. Expanding Eq. (3) in 
Fourier space to first order of the morphological perturbation, we obtain the linearized evolution equation that will 
be used in the following linear analysis: 

dhk{<l, t)/dt = -Tkq^ [lkq^hk{(\, t) + 4(q)] , (4) 

with f/c(q) the first order elastic energy density evaluated at the growing surface position. For the buried interfaces 
between layers z + 1 and i [1 < i < k — \), the morphology hi{x, y) = h + X]q '^i(q) exp(i(7a;X + iqyv) as well 

as its perturbation hi arc also determined by Eqs. (3) and (4) during the growth of layer i, with the parameters of 
layer k replaced by those of layer i; then the morphology at the top of the ith layer remains unchanged (with the 
evaluation at t = li/vi) when the subsequent layers are deposited. 

The above evolution form of the surface morphology is analogous to that used for single-layer film growth, ^"^"^^'^^'^^ 
but the crucial difference is in the determination of elastic energy density £, which is affected by all the underlying 
layers and the substrate, which are made up of different materials, in the inhomogeneous multilayer systems. To 
obtain the elastic free energy in linear elasticity theory, the stress tensor of each layer i (1 < i < fc or substrate i = 
with eo = 0) 



-Upp,iOal3 + Ua0,i - r— CiOa/? 



(5) 



(a, f3 = X, y, z) should be evaluated. Here /ij = Ei/2{1 + fj) is the shear modulus of layer i, with E^ and Vi the 
Young's modulus and Poisson ratio, and Wa/3,j = {9aU/3,i + dfjUa^i)/2 is the linear elastic strain tensor governed by the 
elastic displacement field Ua,i in layer i. 

The ideal way to determine the elastic state of the system is to directly solve the mechanical equilibrium equation 

d0(Jal3,i = 0, (6) 
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and then obtain the elastic fields Ua^i for all the coupling layers and substrate (0 < i < fc) without any approximations. 
Here we have assumed that in this nonequilibrium system of growing multilayer film, the mechanical equilibrium (6) 
occurs instantaneously compared to the film morphological evolution process. The corresponding boundary conditions 
can be obtained by the consideration of the negligible external pressure on top growing surface, leading to 



cra0,kni3,k = at z = hk{x,y,t), 
and the coherence at the buried + interface for < i < — 1: 

<^ai3,i+iri0^i = aai3,ini3^i and = Uq,,, at z = hi{x,y), 



(7) 



(8) 



with Ufs^k and n/3,i the unit vectors normal to the surface of layer k and {i + interface, respectively. Also, in the 
semi-infinite substrate, the decay of strain field far away from the film results in 







for 



-oo. 



(9) 



Therefore, for a fc-layers system with one top surface and k interfaces, wc have to solve 6(fc + 1) equations from 
boundary conditions (7)-(9) to determine the solution of mechanical equilibrium equation (6) for all the layers and 
substrate. Also, due to the nonplanar interfaces between different layers, even if the elastic constants of all the layers 
and substrate are assumed to be identical, the stability results of the system here are not equivalent to those of a 
simplified semi-infinite solid as in the former study of single layer film with flat film-substrate interface. So to 
obtain the exact and direct solution of these equations is very difficult for large fc, even in a linear analysis. 

Here we propose a recursive procedure to directly calculate the strain fields of this multilayer system to first order 
of perturbations. In the framework of isotropic elasticity and the linear analysis, the elastic free energy can be 
determined from the exact solution of the (first-order) elastic fields and then used in dynamical equation (4) for the 
analysis of morphological evolution. 

The basic state of this multilayer structure corresponds to a uniform film with a planar growing surface at constant 
rate Vk and all the planar and coherent interfaces imderneath. Since the basic state of the substrate is stress-free, 
with Ua,o = Ua/3fi = = (a, (3 = X, y, z), and all the surface and interfaces are flat, due to the coherency of 

the film there are no in-plane displacement and strain for all the layers, so that Ux,i = Uy^i = and Ux0,i = Uy^^i = 
(/3 = X, y, z and 1 < z < fc). Also, from the boundary condition (7) on the top surface and the coherency at 
all the flat interfaces, we have aaz,i = with a = x, y, z in each layer i. The displacement in the z direction 
is nonzero: Uz,i = Ui{z — h) + X^i=i ^ih with Ui = ej(l -|- — fi), and the Poisson relaxation in the z 

direction leads to constant strain u^za ~ Ui. The basic-state film is stressed in the lateral direction, with stresses 
^xx,i = ^yy,i = = —"^IJ-iUi for all the layers. Note that in general case, the elastic constants used here (Ei, Vi, and 
/Uj) can be different in different layers and substrate. 

The expansion form of the mechanical equilibrium equation for each layer or substrate can be obtained by substi- 
tuting Eq. (5) into Eq. (6) and expanding in Fourier space with respect to perturbations around the above basic 
state (in the following expressions all the hatted quantities denote the corresponding Fourier perturbations): 



(1 - 2ui){dl - q^) 







iqx 




+ 


iqy 









iiQxUx,i + iqyUy^i + dzUz,i) = 0. 



(10) 



for ith layer (1 < i < fc) or substrate {i — 0), with general solution 
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(11) 



for 1 < i < A:, and 
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for substrate i = due to boundary condition (9), where the parameters Ci = {iqxCtxr 



iqyC 



(12) 



Di = {iqxPx,i + iQyf3y,i + qoLz,i)/{'i - 4i/i), and B = {iqxUxfi + iqyCtyfi + qazfi)/{3 - 4,vo)- These forms of differential 
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equations and general solutions are the same as those for single-layer film.,^° but with different linearized boundary 
conditions. From Eq. (7), we have 

^I3z,k = iqpokhk {P = X, y) and azz,k = 0, (13) 

at top surface z = hk = YIh=i + ^fc*' while at the interface between layers i + 1 and i, i.e., at z = hi = Ij 
(0 < i < fc — 1), boundary condition (8) leads to 

-iq^io-i+i - &i)hi + a/3z,i+i = o-/3z,i {0 = X, y), 

o-zz,i+i = o-zz,i, (14a) 

and 

w/3,i+i = w/j.i {P = x,y), 

{ui+i - Uijhi + Uz,i+i = Uz,i- (14b) 

Note that when i = 0, we have uo = '''o = and ho — C. for the substrate, and condition (9) has already been used to 
obtain the solution form (12). 

All the above mechanical equilibrium equations and boundary conditions are for the /c-layer system, and are also 
fulfilled by the {k — l)-layer system before layer k is deposited, with k replaced by — 1. Let us denote the elastic 
displacement field for this (A; — l)-layer system (with average film thickness h) by Ua^i\k-i with < i < k — 1 

and a = X, y, z, where "|a;-i" represents that the system is of {k — 1) layers. When a new layer k is deposited, the 
elastic state of the whole multilayer is changed and the displacement field is now Ua,i\k (for the fc- layer system with 
average film thickness Y^^Zi h + Vkt). For all the buried {k — 1) layers and substrate, this elastic field is composed of 
two contributions: 

Ua,i\k = K,i\k + K,i\k for 0<i<fc-l. (15) 

Here the first contribution u° ^\k is the old-state field before the deposition of new layer k, and is identical to Ua,i\k-i 
determined by the previous (fc — l)-layer system and assumed to be known. The second contribution u" Jfe represents 
all the change of strain field caused by the new deposition, and is the unknown new-state field. The other unknown 
quantity is the elastic field of the new top layer k: Ua,k\k, and both of them should be calculated from the mechanical 
equilibrium equations and boundary conditions. 

Since the old-state field = Ua.i\k-i obeys the mechanical equilibrium equation (6) (or (10)), from Eq. (15) 

we know that new-state elastic field also follows the same equation, with the same solution forms (11) and (12), 

and so does the unknown elastic field Ua,k\k of the new layer k. For the linearized boundary conditions (13) and (14), 
when we replace Ua,i by formula (15), (Ta/j,, by a-ai3,i\k = <^af3 iU + (0 < « < A; — 1) which can be obtained 

according to Eq. (5), and use the similar boundary conditions for {k — l)-layer system, which are: at 2; = hk-i, 

^f3z,k-i\k-i = iqp^k-ihk-i and ^zz,k-i\k-i = 0, 
and a,t z = hi for {i + l)/i interface (0 < i < fc — 2), 

— iqp{^i+l — ^i)hi + ^0z,i+l\k-l = CFj3z,i\k-l and <Tz2,i+l|fe-l = ^zz,i\k-'L-, 

U0,i+i\k-i = U0,i\k-i and (wj+i - Ui)hi -\- U2,i+i|fe-i = Wz,i|fc-i, 

we can derive the simplified boundary conditions for new-state fields u" Jfe (0 < i < fc — 1) and for the top- layer field 
Ua,k\k- On top surface of the fc-layer system, i.e., z = hk, condition (13) remains unchanged: 

^pz,k\k = iqp^khk {P = x,y) and ^zz,k\k=^, (16) 
while at ^ = hk-i for the k/{k — 1) interface, we have (with (3 = x, y and Ua,k-i\k-i = fe-ilfe for a = x, y, z) 

-iqpakhk-i -\- af3z,k\k = ^}3z,k-i\k and ^zz,k\k = ^zz,k-i\k> (l^a) 
w/3,fe|fe = w/3,fe-i|fe-i + Mjs.fe-iU and 
{uk - Uk-ijhk-i + Uz,k\k = Uz,k-i\k-i + u'^,k-i\k- (17b) 
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At all the other underlying interfaces, the boundary conditions are much simpler: 

Kz,i+i\k = Kz,i\k and ul,i+i\k = K,i\k (18) 

at z = hi with a = x, y, z and < z < — 2. 

Note that the boundary conditions (18) for all the interfaces buried below layer fc — 1 are equivalent to those of 
planar interfaces in the elastic media with displacement fields u^, i\k- Thus, from all the boundary conditions (16)-(18) 
and the mechanical equilibrium equations followed by elastic fields Ua,k\k and w" Jfe, we can obtain an effective elastic 
system for which the elastic state is identical (to first order) to that of the fc-layer system we are interested in. This 
effective system is composed of a growing and undulating surface of top layer k with condition (16), a nonplanar 
interface between layers k and fc — 1 with special boundary conditions (17), as well as the effective planar interfaces 
(with respect to the new-state fields Jfc) between all the other buried layers and between film and substrate. The 
schematic illustration of this procedure is shown in Fig. 1 (a) and (b), where part (a) presents the original fc-layer 
system with nonplanar surface and interfaces and with elastic fields Ua^i\k (0 < z < fc) to be determined, and part 
(b) is the effective system (related to elastic fields Ua,k\k for top layer k and wJJ Ja: for all layers underneath) that we 
obtain here to first order of perturbations. 

This elastic system can be further simplified and then solved exactly when we consider the identical elastic constants 
for all the layers and substrate. The assumption of identical elastic constants is applicable to the system of which 
substrate and all the layers in the film are composed of materials with similar elastic constants, and is used in the 
following calculations so that we could focus on the effects related to surface mobilities and layer thicknesses. Note 
that for a single-layer film, grown on a flat substrate, the elastic behavior of the system is equivalent to that of a 
semi-infinite stressed solid in the case of equal elastic constants. ^"'^^ This can also be obtained by combining the 
elastic solution forms in Eq. (11) for layer 1 and in Eq. (12) for substrate with the flat boundary condition (18), 
which leads to the identical (new-state) elastic fields for the first layer and substrate, i.e., Ua i\k — o\k- Similarly, 
from all the effectively planar boundaries below layer fc — 1 (as shown in Eq. (18) and Fig. 1 (b)), we can obtain the 
same solution for the new-state elastic fields in all the buried layers and substrate: u2,oU = '"a,iU = ■•• — "a,fe-iU- 
This implies that the elastic state (with respect to Jfc) for all these layers and substrate (0 < i < fc — 1) is equivalent 
to that of a single semi-infinite effective medium. 

Consequently, under the approximation of identical elastic constants and in a first order analysis, we can simplify 
the original fc-layer structure to an elastically equivalent system composed of a single layer fc (with elastic field Ua,k\k) 
and an underlying effective semi-infinite elastic medium (with respect to field fc-iU); as illustrated in Fig. 1 (c). 
Now we only have one nonplanar interface at z = hk-i, but with special boundary condition (17) which is different 
from that of the usual single-layer film growth. The undulating free surface at z = hk satisfies the condition (16) 
as usual. The elastic state of this effective system can be solved directly, in a manner similar to that for the case 
of a single-layer film with rough film-substrate interface (but with different boundary condition). Thus, from Eq. 
(15) (and noting the boundary condition (17)), we can obtain the elastic fields of the original fc-layer system Ua,i\k 
{0 < i < k) expressed by that of the previously grown (fc — l)-layer system Ua.j\k~i (0 < j < fc — 1). 

Using the same method, elastic fields Ua.j\k-i can be calculated in terms of Ua,j\k-2 for the (fc — 2)-layer system. 
This procedure is repeated by expressing the elastic fields for z-layer system in terms of those for (i — l)-layer system 
with i = fc, fc — 1, fc — 2, 2, until we reach the single-layer system for which the elastic fields Ua,i\i {i = 0, 1) have 
already been obtained for both the cases of planar^'' and rough^* film-substrate interface. Thus, from this recursive 
procedure we can derive the first-order elastic fields of the fc-layer growing system in the expression of known material 
parameters, and then calculate the elastic free energy for the determination of morphological evolution equation. 



III. MORPHOLOGICAL EVOLUTION FOR GENERAL CASE 



In the linearized equation (4) for surface morphological evolution, the first-order elastic energy density is 5fc(q) = 
^k{uxx\k +Uyy\k) with strain tensors iixx\k and Uyy\k evaluated at top surface z = hk- For the case of identical elastic 
constants {E, ^, and y) in all the layers and the substrate, the method in Sec. II (or the effective system in Fig. 1 
(c)) yields 



4(q) = -E'elqhk + E' 
E 



-qlk 



qhk- 



+ 



1-v 



^ 2(1-1/)' 

{qlk - l)e"«''= (zga;Wx,fe-l|fe-l +'>'qyUy,k-l\k-l)^^h^_^ 



- qhe 9'" (gu^,fc_i|fc_i)^^^^_^ 



(19) 
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where E' = 2E{1 + ly) / (1 — v) and Ik ~ Vkt is defined for the top growing layer k. Also Ua^k-i\k-i{z = /ife-i) is 
obtained from the solution of elastic fields according to the recursive procedure presented in Sec. II, with the result: 



{iqxUx,k-l\k-l +iqyUy,k-l\k-l)z^h^_^ = 2(1-1/) 



quk-ihk- 



k-2 



k-1 

- q y ^ {uj+i -Uj)e hj-quie •=! C 



{iqxUx,k-i\k-i + iqyUy,k-i\k-i - quz,k-l\k-l)^^f^^_^ = quk-ihk-1, 



with Uj = ej{l + i/)/(l — v). 

Substituting Eq. (19) into Eq. (4), we obtain the evolution equation of surface morphological perturbation: 

dhk{q,t)/dt = akhk{q,t)-TkE'ekq''e-i'"'' 

K^.^ ''-1 
^(e,+i-e,-)e *=^+^ hj{q)+eie ({q) 



where 



cTfe = Tk {E'elq^ 



■ ikq 



(20) 



(21) 



corresponds to the (single-layer) perturbation growth rate without the coupling of buried layers and interfaces, and 
the interface perturbations hj and C, are time independent due to the absence of interdiffusion. In Eq. (20), the 
influence of elastic properties and morphologies of all buried layers on the evolution of the top surface profile is shown 
explicitly. Note that due to the exponential decay term exp(— g Y^'lZjj^i k), there is much smaller influence from lower 
buried layers, as usually expected. 

Note that the above linearized dynamical equation for surface morphology is different from that in Ref. 24 by 
Shilkrot et al. for a periodic structure composed of alternating layers of materials A and B (see Eq. (10) there), 
which was also derived to first order of perturbations with the assumption of equal elastic constants within the 
multilayer structure. This is attributed to the different methods in the clastic analysis of the system. Shilkrot et 
al. used the Eshelby procedure^^ by assuming each buried interface as the interface of a misfitting inclusion in a 
semi-infinite homogeneous and isotropic elastic medium with free surface, and considered the contribution from each 
interface to be separate. Here we apply a more direct procedure to determine the elastic state, and find that for this 
highly inhomogeneous system, an underlying interface can only be treated as one between an effective semi-infinite 
homogeneous medium (with respect to a new-state elastic field u^) and another elastic medium (of top layer) with 
finite thickness and free top surface, under the condition of identical elastic constants and linear analysis (see Fig. 1 
(c)). The explicit form of our evolution equation is simpler than that of Ref. 24 and so is our solution as well as the 
corresponding recurrence relation as shown in Sec. IV for periodic A/B system. 

Using the initial condition hk{q,t = 0) = hk-i{q), we can easily solve the evolution equation (20), with the result 
for the top surface morphology: 



hk{q,t) = 



qvk 



+ c^ky 



[{qVk + TkiE'ekCk-iq^ - Ikq")) e'^"* + TkE'ek{ek - ek-i)q^ e-'^^"'] h-i 



+ (e-«'"=* - e"'-') TkE'qhk 



■fe-2 

.i=i 



9 E h- 



fc-1 -I \ 

■'£."<]}. 



(22) 



which is one of the central results of this paper. In Eq. (22), the morphologies of buried interfaces {hj, j — 1, 
A; — 1) remain unchanged during the evolution process of top layer, as a result of negligible bulk diffusion of the film, 
and can be determined from the similar calculation for lower j-layer system. In the case of nonzero misfits for all the 
deposited layers, (i.e., Cj for all i = 1, 2, /c), Eq. (22) can be reduced recursively through the expression of the 
surface or interface perturbation in terms of the perturbations of two lower successive interfaces, leading to the form 



hk{q,t) 
hk-iiq) 



fc-1 

3=3 



h2{q) 
hi{q) 



(23) 
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In Eq. (23), Jjk and Lj (0 < j < fc - 1) are 2 x 2 matrices: 



where 



ifc,i(g, i) = \ { [qvk + Tk {E'ckek-iq^ - jkQ^)] e'^'=* + TkE'ek{ek - ek-i)q^e-'""'*} 

qvk+CTk 

«Ufe + c^fe rfc_iefc_i e-?'"-! - e<^^-i''=-i/'"=-i ^ ' ^ ' 

and Lj_i{q) and Lj,2{q) are obtained by replacing A: with j as well as f with Zj/wj in _L/j,i((/, f) and Lk.2{q,t), respectively. 
Results for first and second layer perturbations hi and h2 are presented in the Appendix, for both rough (C ^ 0) 
and planar {C, = 0) film-substrate interfaces. Therefore, from Eqs. (23) (25) and the results in the Appendix, we 
obtain the explicit expression for first-order surface morphological perturbation hk (<?, t) in terms of known material 
parameters, which can be directly used for stability analysis as shown in the following sections. 

Note that the above solution results apply to the general case for which growth parameters (including the layer 
thicknesses, misfits, surface mobilities, surface tensions, and deposition rate) of different layers can be different, and 
then apply to both the non-periodic and periodic multilayer structures. Also wc have assumed the identical elastic 
constants in the whole multilayer film, and the nonzero misfit strains for all the layers. In the case of zero misfit = 
for some spacer layers i (corresponding to the strained/spacer multilayer structures), the solution forms are simpler 
and will be discussed in Sec. V. 

IV. TENSILE/COMPRESSIVE MULTILAYERS: PERIODIC A/B SYSTEM 

Here we apply the results for morphological evolution, derived in the last section, to the periodic tensile/compressive 
multilayer structures that have been investigated in many experiments, such as InAs/AlAs on InP(OOl),^"^ 
GaInP(InGaAs)/InAsP on InP(OOl),®"* etc. This type of multilayer film consists of two kinds of alternating lay- 
ers with materials A and B. Then we have misfit strain eA, layer thickness Ia, surface mobility Ta, surface tension 
7A, and deposition rate va for odd-number strained layers, and the corresponding parameters es, Ib, ^b, 7b, and vb 
for even-number layers, with nonzero misfits for all the layers: ^ and es 7^ 0. 

A. Formulation of morphological evolution 

For this periodic A/B system, the solution for the top surface morphological perturbation can be directly obtained 
from Eqs. (23)-(25). For a fc-layer structure, when k is odd implying that the top layer k is of material A, we have 

with 2x2 matrices: 



while for even number of deposited layers, the top layer k is of material B and the solution form is 

'f'--)^-L«L.(L.L.) (*>4). (28) 



Here, the solutions (26) and (28) are expressed in dimensionless form, and we have used the rescaling 

q* = ql% T = t/T\, and h* = hi/l\ (0 < i < fc), (29) 
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with the characteristic length and time scale respectively given as: 
and 



(30) 



(31) 



These are determined according to the (single-layer film) growth of layer 1 (assumed of material A). We also introduce 
the characteristic velocity 



■7,0 



(32) 



as well as the following transformations 



'i^*a = '"aIv\, v*g=VBlv\, 

^*B=^bI^A, i*B = ^BltA, l*B=lBhA- 



(33) 



Then the matrix elements in Eq. (27) are rescaled as 



LA2{q*,T) 

LBi{q*,T) 



q ^A^"A 

LA2{q*,T) e-^-B/v'sWs^ 

q*+a*B/v*B 1 



+ (l-e^)g*V''*''^^} 



q*v\ + a\ {T*Blv*B)e*B e 



-q'l*B - e(<^B/''B)'is \ / 



LB2{q*,T) = - 



q* + (7*b/v*b ^ ^ 
q*v\^a\ r 



Vb ^ 



* *4 1 

7b9 



B * 



-— ^ [e-'i^Br _e-yBA 

-g'lA - e'^A^A^A V ) 



with LAi{q*), LA2{q*), LBi{q*), and LB2{q*) obtained by replacing t with Ia/ua or 1*b/v*b in Eq. (34), and 



(34) 



* * 3 *4 

(^A = q -q , 



<^*B = ^B (' 



, 2 »3 

eb q 



ihq 



(35) 



^From Eq. (35), we can have cr^ > or a*B > for a certain wave number q* , and then the expressions in Eq. (34) 
show that |.L^i(2)| ^ 1 or 1^31(2)1 ^ 1 when time r is very large. Consequently, we can obtain from Eqs. (26) and 
(28) that for a certain top surface layer, when the growth time is large enough the amplitude of surface perturbation 
h will grow rapidly with time. Thus, in principle this multilayer system is unstable. 

However, in real experiments the thickness of each layer {I a or Ib) is fixed, including the top layer, then the growth 
of this multilayer system could be effectively stable or unstable?'^ ''^^ That is, when the film surface is flat after the 
growth of large number of strained layers (fc ^ 1), i.e., the initial perturbation decays with further deposition, the 
system is effectively stable; otherwise, if the initial perturbation is amplified and grows in amplitude, the system is 
unstable. For the latter case, the film surface is observed to be very rough, with regular or irregular modulations. 

The stability properties can be calculated from Eqs. (26) and (28), where lt^{q* ,t), "La, and should be replaced 
by hl{q*), Lia, and Lb, respectively, and then the top surface perturbation is mainly determined by the power of 
matrix L = Lb^a, which has two eigenvalues Ai and A2 given by 



if- - - - r- - - -2 - -l^/^l 

M,2{q*) ~ 2\ (^^"^ -^-^2 + LaiLbi) ± y{LA2 + -£'B2 + LaiLbi) — 4L^2-^^B2 } ■ 



(36) 



Assuming that the corresponding two eigenvectors of matrix L are linearly independent, we obtain the explicit 



10 



expression for the top surface morphological perturbation h'^{q*): 



h*k = Lai^-t h2 + LA2 -Lb2^-: 

Ai — A2 \ Ai — A2 




+ 




k- 



- 1 k- 



) 



hi {k odd), 



( 




k k k-2 k~2 



) 



^2 + L1A2LBI 



^1 ^ ~ ^2^ 

Ai — A2 



fc-2 fc-2 



/ij (fc even) 



(37) 



which can be easily proved to be proportional to A*^/^ [A = niax(Ai,A2)] for real eigenvalues, or to |A|*/^ (|A| = 
I All = IA2I) for a complex eigenvalue pair, when the number of total layers A; ^ 00. Thus, the effective stability of 
this multilayer system (defined for k ^ 1) is determined by the magnitude of eigenvalue |A|. When the maximum 
|A| (with respect to wave number q*) is larger than 1, the surface perturbation grows, leading to the instability 
of the multilayer film; otherwise, if |A| < 1 for all the q* modes, the system is stable due to the decay of surface 
perturbations. According to Eqs. (34) and (36), the value of |A|, and then the stability properties of the system, 
depend on 6 dimensionless parameters Tq/vq, e*^, v'^, 1^, and 7^ that are defined in Eq. (33). 

The expression for A remains unchanged with respect to the interchange of A and B layers, since the eigenvalues of 
matrices L = LbLa and L' = L/iLg, with liAiB) defined in Eq. (27), are the same. Thus, the stability properties of 
multilayer structure are determined by |A| and are identical in the k 00 limit for the A/B and B/A stacking forms. 

Note that the stability analysis presented here is for very large fc 1, while for the morphology of each layer 
i (i = 1, fc), or for the phase^'' between two successive strained layers, which depends on the sign of h*/h*_i, 
the evolution Eqs. (26) and (28) [with the explicit expressions similar to Eq. (37)] should be used. Although the 
stability properties for k ^ 1 are independent of the type of the first layer or the; interchange of A and B layers in the 
structure, the detailed morphology of each layer i may be influenced by whether the first deposited layer is tensile or 
compressive, as seen from the expressions (26), (28), and (37). 



In this subsection, we present the results of stability diagrams for various (rescaled) material parameters and growth 
conditions. The stability properties are determined through the maximum magnitude of eigenvalue A, with the use of 
the above analytic expressions and discussions. One of the factors that wc focus on is the surface mobilities of strained 
layers. During the growth process, different mobilities of A and B layers results in different interface morphologies 
which are frozen in due to the deposition of upper layers, but which can influence the system's strain field and the 
evolution of top surface profile. 

A common feature of the stability diagrams presented in this subsection is that the multilayer system can be 
effectively stabilized by large enough deposition rate (u^ and w^), as also found in experiments.^'* This is similar to 
single-layer film growth:^°"^^'^'''^* for large enough deposition rate, the morphological perturbation of layer surface 
does not have enough time to develop through surface diffusion, since it is buried and frozen by fast subsequent 
materials deposition. 



Stability diagrams in the /^-wj space are shown in Fig. 2 for different relative surface mobilities = Fb/F^ = 
0.001, 0.01, 1, and 100 (other parameters are e|j = —1, l\ = Ig, = Vg, and 7^ = 1). Important features arc: 
larger deposition rate imply greater stability, very thick or very thin layers imply instability and for large I a and I b , 
multilayer system loses coherency; larger relative mobility F^ leads to more instability (note the larger unit-scale for 
v\ in the inset). 

To understand the effect due to different mobilities, we consider that the increase of F^ = F^/F^ corresponds 
to keeping the parameters of material A (including Ta, ^a, and 7^) unchanged, as seen in the rescaling Eq. (33), 
and to the increase of surface mobility Fs of layer B. Consequently, for larger F^, the morphological instability is 
enhanced due to a faster diffusion on surface of layer B (driven by nonuniform strain energy distribution^^'^"), and 
then larger deposition rates v*^ and v*g are needed to suppress it. In multilayer experiments, different mobilities result 
in different morphological profiles of A and B layers.^ Layers with high mobility exhibit the morphology of ripples 
with high peaks and deep valleys, while for layers with low mobility, flatter profile is observed. Therefore, increasing 
the surface mobility of either type of layer leads to rougher surface profile, which is more difflcult to be smoothed out 
during subsequent growth. 



B. Stability results 



1. Effect of different surface mobilities in strained layers 
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Our analytic expressions in last subsection can also explain this effect. Equations (36) and (34) show that the value 
of |A|, which determines the effective stability of the system, depends on parameter Tg/v"^. Thus, when r|j is larger, 
deposition rate should be correspondingly increased in order to keep the same value of r^/t;^, or to retain the 
system within the region of stability. 

2. Interplay between mobility and global misfit of multilayer 

In the tensile/compressive multilaycir structures, one important parameter is the global strain of the whole multilayer 
film with respect to substrate, which varies with the relative layer thickness Ib/^a {— Ib/Ia) and relative misfit 
(= es/f-A)- A rjualitative measure of the global strain is provided by Ia^a+Ib^b- It is interesting to study the stability 
properties of systems with different global misfits, in particular the asymmetry for global tensile and compressive 
multilayers and its relationship to different surface mobilities of composite layers. 

The effect of relative thickness Ib/1*a (with fixed relative misfit = —1) is depicted in Figs. 3 and 4. In Fig. 3 (a) 
we plot the stability diagrams of l\ versus v\ for different values of Ig/l'X, with small relative mobility Tg = 0.01. 
When increases from 1, the stability boundary shifts to the lower value of l\\ but when 1*b/1\ < 1, the system 

is much more unstable, that is, very high deposition rate v\ (orders of magnitude larger than the scale shown in the 
figure) is needed to stabilize the multilayer morphology. 

Stability diagrams as a function of 1*b/1*a and v\ (where the A/B period thickness l\ + Ig is held constant as in 
experiments) can illustrate this asymmetry more clearly, as shown in Fig. 4 (a) and (b). When Tg < 1 (> 1), i.e., 
Fb < (>) Fa, the systems with Ib/Ia > 1 (< 1) are much more stable, corresponding to much smaller deposition 
rate v\ (= Vg here). When r|j = 1 (F^ = F^), Fig. 4 (b) shows that the stability diagram is much more symmetric. 
Also, for different relative mobilities T*b [e.g., 0.01 and 0.001 as in Fig. 4 (a)], the system can display cither a greater 
instability or a greater stability as the layer period thickness Ia + ^b increased. This can be understood through the 
stability diagram of l\ versus v\ in Fig. 2, where within different range of values of l\, one can see that the stability 
can cither increase or decrease with layer thickness This stability property with respect to period thickness is 
richer compared to the observation in InAs/AlAs experiments, where the modulation amplitude of the films was found 
to decrease with the decrease of period thickness.^ 

Figures 3 (b) and 5 exhibit the effect of relative misfit e^, with fixed 1*b/Ia ~ "^^^ asymmetric property, seen in 
Fig. 4, also appears in the stability diagrams of versus for different values of eg [Fig. 3 (b)]. For F^ = 0.01, 
when |e^| > 1 the instability is enhanced by the increase of the value of misfit, with the shrinking of stability region; 
however, when |e^| < 1 the v'^ values of stability boundary are much larger than the scale shown here, corresponding 
to a very unstable system. We can see the asymmetry more clearly in the e|j-w^ plot of Fig. 5, where Ia^Ib~ 0-5 
and when F|j < 1 (> 1), the systems with larger (smaller) misfit of layer B [i.e., |e|j| > 1 (< 1)] are much more stable. 

Consequently, these stability diagrams in Figs. 3-5 show that when the sign of the multilayer global strain is the 
same as that of the layers with smaller surface atomic mobility, the system is more stable; otherwise, when the layers 
with larger mobility have relatively larger misfit strain or layer thickness, the modulations on the surface of these 
layers are enhanced due to faster surface diffusion and arc hard to be suppressed by the subsequent deposition of 
layers with smaller mobility. Note that this result of stability asymmetry, which looks different for F^ > 1 and < 1, 
is in fact symmetric with respect to the exchange of A and B layers. E.g., when layer A is compressive, layer B is 
tensile, and F^ = Tb/Fa < 1, we obtain from above analysis that global tensile multilayer film is more stable. If A 
and B are interchanged, layer A' (of material B) is tensile and layer B' (of material A) is compressive, with F|j, > 1; 
then according to above result, the global strain of the more stable system should be the sign of layer A', which is 
also tensile. 

In recent experiments of InAs/AlAs short-period superlatticcs grown on InP(OOl),^"''' the system with large tensile 
global strain was observed to be stable. This is consistent with our theoretical results, since AlAs is tensile with 
respect to substrate InP, InAs is compressive, and importantly, Al has much smaller surface atomic mobility, leading 
to the possibility of stabilization for global tensile multilayer (corresponding to thicker AlAs layers) due to the above 
analysis. For the global compressive films (corresponding to thicker InAs layers), our calculations yield much larger 
instabilities, in agreement with the most recent experimental observation^ that although the multilayer structure 
under global compression exhibits weak lateral modulations as a whole, the system is in fact morphologically unstable 
due to a high degree of surface roughness and the reduced strength of lateral modulations can be attributed to irregular 
surface patterns. 

To further compare with the experiments of InAs/ AlAs, stability diagrams of Ib/Ia versus (Fig. 4) should 
be used. Here layer A corresponds to InAs and layer B to AlAs. As shown in the figure, the stability boundary 
depends on the parameter values, in particular the relative mobility F|j (which also depends on growth temperature 
as discussed below: when F|j < 1 that applies here, higher temperature T leads to larger r|j) and period thickness 
l\ + Ib- When F^ is very small (e.g., = 10~^) and layer thickness l\ + Ig is thin enough (e.g., = 0.3), our results in 
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Fig. 4 show that films with larger tensile global strain (corresponding to larger Ig/l'X) are more stable. However, when 
layer thickness + Ig increases (e.g., = 0.9), we have opposite result for the multilayer under global tension; higher 
global strain leads to more instability. When relative mobility is larger (e.g., T*g = 10^^), difTcrcnt behavior can be 
obtained. Opposite to the case of r|j = 10""^, thinner A/B layers (Z^ + = 0.3 in Fig. 4) result in more unstable 
global tensile film. Interestingly, for larger thickness 1% + ^% — O-^? from Fig. 4 we can get a maximum boundary value 
of (maximum instability) in the tensile region when Ib/Ia is somewhat larger than 1 (i.e., the whole multilayer is 
somewhat tensile), which may correspond to the experimental observation of the maximum modulation for slightly 
global tensile InAs/AlAs film.^'^ Also there is a minimum boundary value for + = 0.9, showing the maximum 
stability at some higher global strain. 

Although our theoretical results qualitatively agrees with some experimental observations, for the quantitative 
comparison more realistic factors should be considered in our model. One important factor is the "wetting effect" 
especially for the multilayer structures like short-period supcrlattices where the thickness of each strained layer is very 
thin and then the film-substrate and layer-layer interface energies have to be carefully included in section II. This 
effect is absent in our calculation but applicable to the InAs/AlAs experimental system. 

3. Strain-balanced condition 

In some experiments for multilayer film growth, the strains caused by tensile and compressive layers are balanced. 
That is, the system fulfills the strain-balanced condition Ia^a + Ib^-b = 0, or equivalently, 1*b/1*a = ~^h*B-: which 
makes the global misfit of multilayer film related to substrate equal to zero. One would expect that the corresponding 
stability diagrams are less asymmetric compared to the globally strained system studied above. This is seen in Fig. 6 
for diagrams of 1%/Ia versus (with fixed Ia^^b)- However, the asymmetry is still obvious, and from these stability 
diagrams of zero global strain we obtain an interesting property that for all values of relative mobilities F^ = Tb/Ta, 
the multilayer structures are more stable, that is, can be stabilized by lower rescaled deposition rates v"^ and v%, 
when the relative misfit |e|j| = |eB|/|eA| < 1 (or equivalently, Ib/Ia > !)• Figure 6 should be compared to Fig. 5 (c) 
and (d), since the later does not have a strain-balanced condition. 

Note that this result is obtained for fixed material parameters F^, f-A and 'ja (especially for rescaled deposition rate 
^A(B) ~ ''^a(b)/va = VA{B)l'A/(^AE'^e^) that determines the degree of stability here), and one can obtain the same 
result with respect to fixed B parameters by interchanging A and B in the calculation (see also similar diagrams for 
Fg < 1 and > 1). Thus, when we already have one layer material 1 (with misfit ei), to make the multilayer system 
more stable, the other layer material 2 should be selected to have smaller misfit in relation to the substrate (i.e., 
|e2| < which implies a larger layer thickness I2 > h due to the strain-balanced condition. 

Consequently, with the condition of strain balance, the multilayer structure with asymmetric A and B layers 
(consisting of alternately thin strained layers and thick layers with less misfit of opposite sign), are more stable 
compared to that with symmetric properties ca = — es and Ia = Ib, and more the asymmetry, more the stability. 
Also, this result is independent of the relative surface mobility of two types of layers. This theoretical prediction, which 
is different from the case of nonzero global strain studied above, should be related to different multilayer structures 
with different types of composite layers, and is yet to be verified by experiments. 

In fact, this result can lead to the well-known fact of the heteroepitaxial growth: the system is more stable for less 
strained composite layers. When we fix material A and select material B, the system is more stable for the material 
with smaller misfit cb (and then larger thickness Ib) according to our result; then, if we fix material B and select 
other type of layer A, the material with smaller eA corresponds to more stable system. The repeat of this selection 
procedure leads to an obvious result that the most stable multilayer system is the structure in which both misfits eA 
and Cb are as small as possible. 

4- Role of surface tension and growth temperature 

Our results for stability properties are not very sensitively dependent on relative surface tension 7^ =7^/7^ (One 
of the major reasons may be that wetting effect^'' is not considered here). See diagrams of versus in Fig. 7 (a), 
with F^ = 0.01 and different 7^ = 0.1, 1, and 10. For relative tension 7^ < 1, the stable region is smaller; but the 
stability diagrams do not change much if we compare the result of 7^ = 10 to that of 0.1 (with the difference of two 

orders of magnitude for "/g). 

We also calculate the diagrams as a function of 1*b/1*a and v*^ in the strain-balanced condition, for different relative 
surface tension 7^, as shown in Fig. 7 (b). The results depend on the relative mobility F^ and period thickness 
l*A-\-l*B- For l*A + l*B = 0.9 and F^ = 1, the inset of Fig. 7 (b) shows that the stability of the system increases with the 
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increase of 7^, while for smaller = 0.01, the results are more complicated: for large Ib/Ia, 7s = 1 (i-e., "fA = 7b) 
corresponds to the most stable case, but for small Ib/Ia, larger "/g favours the stability. 

Growth temperature is an important factor in real experiments. Our results in Sec. IV A show that the stability 
properties are affected by growth temperature T via mobilities Ta and Tb, since the rescaled parameter changes 
exponentially with T: oc exp[— — i?^)/fcBT] with Ea and Eb the activation energies of A and B materials. If 
Eb > Ea, as temperature T increases, value of Tg increases or remains unchanged; while for Eb < Ea, r|j decreases 
with the increment of T. Also, the dimensionless deposition rates Va and v'^ decrease with temperature due to 
^A(B) = VA(B)lA/{'^AE'h\) [see Eqs. (32) and (33)], Fa oc T-^ exp{-EA/kBT), and then v*^^j^^ cx T exp{EA/kBT). 
Thus, when temperature T increases, for the case of Eb > Ea (i.e., r|j < 1) larger (or fixed) value of Tg and smaller 
deposition rate v*^^g^ render the system more unstable, according to the stability diagrams in e.g.. Fig. 2. Similar 
result can be obtained for Eb < Ea (i.e., T*g > 1). As discussed above, when T increases, value of ^^(5) gets 
smaller, leading to more unstable system, but F|j also decreases, favoring the stability of the film. Note that due to 
< Ea — Eb < Ea, the reduction of F^ is much slower than that of v'^f^g^ as T is increased. Hence, we still have the 
result that higher growth temperature enhances the instability of this multilayer structure, as expected. 

C. Discussion 

For the alternating tensile/compressive multilayer, when a layer ?', which is tensile (compressive) with respect to 
the substrate, is grown with the presence of undulating surface, local strain energy concentration is created at the 
surface troughs with more stress relaxation at the peaks. When the next compressive (tensile) layer i + 1 containing 
larger (smaller) atoms is deposited on it, the difference of strain energy density along the underlying surface drives 
the deposited atoms to the trough regions. This corresponds to the coupling of strain fields for two successive upper 
layers, while in reality the strain coupling is more complicated due to the influence of lower buried layers, which has 
been shown to decay exponentially with the distance from the growing surface [see Eq. (20)]. Nevertheless, some of 
our results in Sec. IV B can be explained by this simple two-layer coupling picture. When the layer thickness is too 
small, that is, the deposited material is not enough to smooth out the underlying troughs, the subsequent deposited 
layer i + 2 will deepen the corrugations since it is of the same material as layer i. Thus, the surface perturbation is 
amplified by subsequent depositions of different layers, rendering the system unstable. On the other hand, when the 
layer thickness is too large, we have inverse morphology, that is, peak (valley) forms and continues to grow above 
the valley (peak) of underlying layer i, and then the surface morphological modulations can not decay with the layer 
stacking. Thus, for intermediate range of layer thickness, the multilayer is more stable, i.e., lower deposition rate 
is needed to reduce and suppress the surface undulations during each layer growth. Here, surface mobility plays 
an important role. Smaller mobility leads to slower stress-driven surface diffusion process, and then the surface 
morphology has less dependence on the corrugations of underlying layer. Consequently, the stability results in Fig. 
2 can be reproduced, that is, the system is more unstable for too small or too large a layer thickness, higher relative 
mobility, and smaller deposition rate. 

More interesting results are found for the multilayer films with global misfit with respect to the substrate. When 
the layers with smaller surface mobility have larger thickness or larger magnitude of misfit, that is, the sign of the 
multilayer global strain is the same as that of this type of layers, the system is more stable. This can also be seen in 
the above approximate picture of two-layer coupling. For example, for the growth of layer with small mobility, the 
surface evolution does not depend much on the undulations of the layer underneath, and then with large thickness the 
surface profile is fiatter than that of underlying layer (as observed in experiments); when the next layer with larger 
mobility is deposited, the surface corrugation develops but is limited due to the smaller layer thickness. Therefore, 
the modulations are easier to suppress during the layer stacking procedure. This asymmetric property is consistent 
with the recent experimental results of InAs/AlAs on InP(OOl), where the films imder high global tension have 
been found to be stable^"^ while those under global compression are much more unstable.^ Note that the formation 
of lateral modulations observed in the whole multilayer film is due to the regular surface modulated pattern; for 
irregular surface pattern weak lateral modulations are observed, but the surface could be very rough, corresponding 
to unstable system.^ 

The other property we obtain for the tensile/compressive multilayer is related to the films with strain-balanced 
condition, i.e., zero global strain. We find that with this constraint, the multilayer structure is more stable if one 
type of strained layers is thicker and then has less misfit of opposite sign compared to the other type. To see this in 
cxpc;rimc!nts. one shoiild compare multilayer films composed of different materials (such as A/B and A/B' structures 
both with zero global strain, but Ib ^ I'g and cb 7^ cb')- discussed in Sec. IVB3, this asymmetry property in fact 
leads to the well-known result that smaller misfit of each composite layer results in more stable system. 
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V. STRAINED/SPACER MULTILAYER STRUCTURES 

The other kind of multilayer that we are interested in is the structure with strained/spacer stacked layers. In this 
system, the strained layers (cither tensile or compressive with respect to the substrate) and spacer layers are alternately 
deposited on a substrate of material B, such as lnAs/GaAs(001),"-" Ge/Si(001)/'*-i'^ and SiGe/Si(001).3'iO'is The 
strained layers are odd-numbered {i = 1, 3, ...), with nonzero misfit strains e, ^ 0, while the even-numbered spacer 
layers (j = 2, 4, ...) are of substrate material B, with zero misfit = = 0. Here, we consider the case in which 
all the spacer layers are thick enough, such that the interface between a spacer layer j and a strained layer (j ' + 1) is 
atomically flat [hj =0 for even j) since the growth front of previously deposited strained layer is smoothed out. [If 
in the solution (22), one sets Cfe = e_B = 0, then is negative and for thick enough spacer layers, Eq. (22) leads to a 
flat top surface of the spacer layer. This is also seen in many experiments.] 

The experimental interest for this multilayer structure is different from that of tensile/compressive multilayer due 
to the presence of spacer layers. This system is morphologically unstable (as verified below) , and 3D islands develop 
for each strained layer with the self-organized formation of multisheet arrays of islands along the growth direction. It 
has been observed that both size and regularity of islands increase with the increment of deposition layers, and then 
saturate for large enough layer numbers. In the linear analysis presented here, this property of the island size can 
be explained through the study of kinetic critical thickness for the onset of 3D instability, with details shown in Sec. 
VB. 



A. Elastic analysis and morphological evolution 

The recursive procedure for elastic analysis presented in Sec. II can be applied to the strained/spacer multilayer 
structures. The corresponding results are simpler due to the zero misfit of spacer layers and the flat strained/spacer 
interfaces (i.e., Cj = 0, hj — for even j), and are different for the odd-numbered (strained) and even-numbered 
(spacer) layers. Consider a /c-layer system with top layer k strained with respect to the substrate (i.e., with odd 
number k). Thus, the interfaces between odd- and even-numbered layers (that is, — 1) interfaces with odd i and 
3 < j < /j) in Fig. 1 (a) are planar. The effective boundary conditions at the top surface z = hk and the interfaces 
z = hi (0 < i < fc — 2) below layer (fc — 1) are the same as Eqs. (16) and (18), while at planar k/{k — 1) interface the 
special boundary condition (17) should be revised as {a = x, y, z) 

(^az.k\k = ^az,k-l\k ^ud Wa./c |fe = "a,*;-! U-1 + "a^fe-l U' (38) 

Similar to the discussion in Sec. II, we obtain an elastically effective system composed of an undulating top free 
surface and a planar interface between top single layer k and an effective semi-infinite elastic medium of field u^k-\\k 
with boundary condition (38), to first order of perturbations and with the assumption of identical elastic constants. 
To determine the elastic fields recursively, we need to consider the system with smaller number of layers. For a m-layer 
system with 1 < m < fc, the above results can be directly applied with the replacement of fc by m if layer number m 
is odd. However, for even number m, the results are somewhat different due to the fiat top surface of even-numbered 
spacer layer m. The corresponding effective system is in fact composed of top layer m, with a planar free surface, and 
underlying effective semi-infinite medium of elastic field m-i\mi with a nonplanar buried interface. The boundary 
conditions are different from those of Eqs. (16) and (38): At planar top free surface z = hm, 

while at rough interface z = hm-i, 



'az,m\m — '-' az.m—W'mi 



(40a) 



Ul3,m\m — 'W/3,m— l|m— 1 + ^^,m-ll»" &'^<^ 
-Um-lhm-1 + tiz.mlm = Uz,m-l\m-l + (40b) 

(Here a = x, y, z and /? = x. y.) 

According to the procedure in sections II and III, we can determine the elastic energy density and then derive the 
evolution equation for surface morphological perturbation of top strained layer fc: 

dhk{q,t)/dt = akhk{q,t) + E'rkekq^e-'""'' 

k-1 

-a T U . 

(41) 



E 



k-i k-i 

-q E h . -q E h " 

e e2j-ih2j-i{q) - e •=! eiC(g) 



15 



which in fact can be obtained from Eq. (20) of general case by considering ej = and hj = for spacer layer j. 
For simplicity, here we assume the planar film-substrate interface, i.e., C = 0) ^^'^ the initial perturbation on the top 
surfaces of substrate and of all the spacer layers, from which the strained layers start growing, to be hP{q), which can 
arise from the background roughness due to the noise of deposition flux even if the starting surface is perfectly flat. 
Then the solution of evolution equation (41) is expressed as 



hk{q,t) = Ak{q,t)h\q), 



(42) 



where Ak{q,t) is the morphological perturbation amplitude for top surface of strained layer k. For k = 1, Ai{q,t) 
exp((Tii) according to Eq. (A.l) and the result of single-layer film growth; for fc = 3, we have 



qv3 + a3 ^ ' 



(43) 



and for A; > 5, 



A^{q,t) = Lk,i{q,t) + Lk,2{q,t){Lk-2,i{q) 

fc-5 
2 



1=2 



n -^2^-^1,2(9)1 L2i+i,i{q) 

j=i+i 



• fc-3 

n^2i+i,2(9)| Aaiq) 
i=2 



}■ 



(44) 



where 



Lk,i{q,t) = e' 



Lk,2{q,t) = 



^k^k qVk-2 + O-k-2 



^k-2ek-2 qvk+cTk e'^fc-2'fe-2M-2 - e-''^-^ 



a'^k-2lk-2/vk-2—q{lk-l+lk-2) 



-glk- 



^k-2ek-2 qvk + (7fe e'^fc-^'^-^M-^ - e-^'^-^ 
{qVk-2 - Vk-2lk-2q^) e-^^"-' +Tk-2E'el_2q''e'"'-'^''-'^'"'-' 



(45) 



and ^3(9), ^2^-1-1,1(9), and Z2j+i,2(9) [2 < j < (fc — 3)/2] are obtained by replacing k with 2j + l and t with l2j+i/v2j+i 
in A^{q,t), Lk^i{q,t), and Lk,2{q,t), respectively. Formula (45) shows that when the underlying spacer layer is very 
thick, i.e., Ik-i ^ 1, we have Lk^i(q,t) exp(c7fet) and Lk^2{q,t) — * 0. Then the single-layer solution similar to Eq. 
(A.l): hk{q,t) = exTp{<Tkt)hP{q) is recovered for this multilayer structure with any value of the layer number k, as 
expected. 

Consequently, we have derived the explicit solution [Eqs. (42)-(45)] for the surface morphological perturbation, in 
terms of material and growth parameters for all the deposited layers (1 < i < k): F,, e^, 7t, k, Vi, and the elastic 
constant E' . Also, the above expressions can be directly used to determine the early evolution and stability properties 
of stacked strained/spacer structures, for which the composite materials and growth conditions of different layers 
could be different (non-periodic structure) or the same (periodic structure). 

B. Results: kinetic critical thickness 



Here we apply the general results derived above to periodic strained/spacer structure studied in many 



experiments, 



for which, all the strained layers (odd-numbered) consist of the same material A, with the same 



deposition rate va and average thickness I a, and the thickness of each spacer layer (even- numbered) is fixed as Ib- 
Thus, from Eq. (44) the perturbation amplitude Ak (for odd k and fc > 5) is simplified as 



Akiq,t) = Li{q,t) + L2{q,t) 



Li{q) 



1-L2(g) 2 
1 - L2{q) 



+ L2{q) — A3{q) 



(46) 



where Ln2){q,t) corresponds to -/jfc,i(2) (<?, i) hi formula (45) in which one has set ai = a a, Fj = F^, = ca, 7i = 7A) 
li ~ Ia, and Vi — va {i = k, k — 2), and Li(2){q) = -^1(2) t = Ia/va), so that all of them are independent of layer 
number. For smaller fc, the amplitudes are Ax{q,t) = exp(CTyii) for single-layer film, and 



TAE'e\q^ 

qVA + CFA 



As{q,t) = e-'^' + 



qlB ^^a-At _ ^-qvAt 



(47) 
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as obtained directly from Eq. (43). 

^From Eqs. (46) and (47) wo find that when the growth time is large enough, i.e., t ^ I, the perturbation amplitude 
Ak{q,t) is proportional to exp(cr^t), with positive perturbation growth rate for small wavenumber q as in the case 
of single-layer film. In Fig. 8, we show the maximum value of perturbation amplitude (with respect to wavenumber 
q) A'^^^{t) as a function of growth time t for different layer numbers, and note the positive growth rate. Thus, 
similar to the tensile/compressive structure discussed in Sec. IV, essentially this strained/spacer multilayer system is 
morphologically unstable. However, different from the tensile/compressive films, the system here cannot be effectively 
stabilized (with the decay of initial morphological perturbation) in certain range of growth parameters. This is due 
to the fact that the surface perturbation amplitude always increases with layer number k, i.e., Ak{q,t) > Ak-2{q,t) > 
... > ^3(9, t) > Ai{q, t) > for certain wavenumber q and growth time t (= Ia/va)-, which can be proved analytically 
from Eqs. (46) and (47) (see also numerical result of maximum perturbation amplitude shown in Fig. 8), and can be 
attributed to the coupling of interface perturbations of all the buried layers exhibited in Eq. (41). Also, if the spacer 
layer thickness Ib is large enough that L2{q) < 1 for any mode q, we find the asymptotic expression for amplitude 
Ak, as fc — > 00, as 

Ak{q,t) ^ Li{q,t)+L2{q,t)Li{q)/[l - Uiq)]. (48) 

Nevertheless, during the growth process the surface instability of top strained layer can be kinetically stabilized 
up to a critical thickness, as in the single-layer film growth. ^° For growing layer below this kinetic critical thickness, 
that is, at early time of layer deposition, the surface perturbation does not develop fast enough compared to the 
deposition of materials flux, and then is not observable. The perturbation becomes apparent relative to the background 
disturbance only when the layer is sufficiently thick and then the initial undulation /i° for the growth of top strained 
layer is substantially amplified. As in the study of single-layer film,^^ we use ln{h/h°) 1 as the amplification limit, 
or cquivalently, we use 

ln[Ar"(i -tc)] = 1 (49) 

with tc = hc\k/vA, as a measure of the kinetic critical thickness hc\k for a fc-layer system (see Fig. 8). In Eq. 
(49), we have assumed that the surface morphological property is dominated by the mode of maximum perturbation 
amplitude {A'^'^^), that is, by the wavelength of the strongest instability which, in single- layer case, may correspond 
to the wavelength of the surface roughness pattern observed in real experiments.^" ''^ 

The rescaled kinetic critical thickness hc\k/hc\k=i (with /ic|fe=i the result for single-layer growth) is plotted as a 
function of layer number k in Fig. 9, according to Eqs. (46) and (49). We choose the material parameters of Ge/Si(001) 
multilayer growth,^^'^^ with strained layer thickness Ia = ^ ML, deposition rate va = ^ ML/min, temperature T = 550 
°C, and different spacer layer thickness Ib = 8 nm, 8.5 nm, and 9.5 nm. (Typical Ge layers are rather thin and 
usually expressed in units of monolayers, while the Si spacer layers are relatively thick and appropriately expressed in 
nanometers, as in experiments. ^''"^^ We use the conversion 1 ML = 0.1457 nm for Ge monolayers^'"' in our calculations.) 
For surface mobility Tk , the definitive measurement of its absolute value over a range of temperature is lacking, and 
we determine it from the expression of single-layer result: hc\k=i = 256i'A7i/(27rAe^£"''), as well as the observed 
value of hc\k=i, corresponding to the experimentally measured 2D to 3D transition thickness, in Ge/Si(001) growth. 
As illustrated in Fig. 9, the kinetic critical thickness hc\k first decreases monotonically and rapidly as the number k 
of deposited layers increases, and then slowly until it converges to a saturation value at large enough number k. This 
convergence effect is due to the asymptotic result of perturbation amplitude Ak [see Eq. (48)] when the spacer layer 
is sufficiently thick and then £2(9) < 1- Also, from the comparison of different symbols (corresponding to different 
spacer layer thicknesses Ib) in Fig. 9, we find that hc\k is reduced more rapidly for thinner spacer layer, approaching 
a lower saturation value at larger number k. This can be attributed to the stronger coupling between top surface and 
the underlying layers and interfaces for smaller Ib, as seen in Eq. (41). 

These results shown in Fig. 9, both for the behavior of hc\k versus k and for the effect of spacer thickness Ib, are 
in good agreement with the recent experimental observations^^"^® of critical thickness for 2D-3D growth transition. 
The properties of kinetic critical thickness presented here can explain some experimental results for islands (quantum 
dots) growth. When the thickness of growing film exceeds hc\k, the morphological instability occurs, exhibiting as 
surface ripples or corrugations which have been observed as the precursor to the coherent islands formation. ■^"■'^^ Thus, 
smaller hc\k not only corresponds to the earlier occurrence of islands during growth, but also results in larger island 
size since it implies more time for the development of 3D islands if the nominal thickness of top strained layer is fixed. 
Consequently, the results in Fig. 9 correspond to the experimental phenomenon of the increase and then saturation 
of islands size in upper strained layers. ^'^'^^'^^ 

To further study the effect of layer thickness, we plot the rescaled critical thickness hc\k/hc\k=i as a function of 
spacer layer thickness Ib in the inset of Fig. 9, for different layer numbers fc = 3, 5, 7, and 11. hc\k increases as the 
spacer thickness Ib increases and the layers number k decreases, in agreement with Fig. 9 and recent Ge/Si(001) 
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experiments.^^ When Ib is large enough, the increase of hc\k becomes much slower, with the value approaching that 
of the single layer growth, as analyzed in Sec. VA. The relationship between the rescaled h^k and the thickness 
Ia of buried strained layers is presented in Fig. 10, which shows the opposite effect of Ia- Larger Ia causes rougher 
morphology of buried interfaces, and hence the enhanced perturbation on surface as a result of the strain fields 
coupling [see also Eq. (41)]. Consequently, for larger I a instability occurs earlier, with smaller value of /ic|fc- 

The behavior of kinetic critical thickness hc\k with respect to misfit strain ca and deposition rate va is shown in 
Fig. 11 (a) and (b). The results of A: = 1 there correspond to those of single-layer growing film with single component, 
for which the kinetic critical thickness is found to follow the e^^ power law and increase linearly with va-^'^ Similar 
properties can also be obtained in our calculations for strained/spacer multilayers with layer number fc > 1 and at 
small misfit [see the range ca < 4% in Fig. 11 (a)] and large deposition rate [see Fig. 11 (b) and its inset]. However, 
small deviations do occur for large ca and very small va, as illustrated in Fig. 11, especially for large k systems. 
Aroimd and beyond the misfit of 4%, Fig. 11 (a) shows that hc\k for fc > 1 decreases faster than e^* as misfit increases, 
with the behavior deviated from the power law; also, as presented in Fig. 11 (b), the linear relationship between hc\k 
and VA is deviated slightly at very small deposition rate. These deviations are more obvious for larger number fc. But 
when k is large enough (fc > 11 for the parameters of Fig. 11), the graphs of hc\k for different fc are almost identical 
in Fig. 11 (a) and (b), due to the convergence effect obtained analytically in Eq. (48) and shown numerically in Fig. 
9 for large fc. 

Note that the kinetic critical thickness hc\k studied here, which is for the onset of morphological instability during 
layer growth, is different from the wetting layer thickness /iwL which corresponds to equilibrium critical thickness^^ 
(postgrown and thermodynamically stable) for the 3D island formation. Usually hc\k > ^WL due to the consumption 
of wetting layer in the formation of 3D islands. E.g., in the molecular-beam epitaxy (MBE) experiment of single layer 
Ge on Si(100),^2 hc\k=i = 3.7 ML at temperature T = 700 °C, while hwh = 3 ML. 

For this stacked strained/spacer multilayer, due to the presence of spacer layers which are of the substrate material, 
the mechanism for morphological instability is more analogous to that of single-layer film growth compared to the 
tensile/compressive structure. But the crucial factor is still the coupling of strain fields in different strained layers. 
Also, from the general expressions in Sec. Ill, it is straightforward to obtain the results of morphological perturbation 
for the case of nonplanar and rough spacer layer B, that is, = but hi ^ for even number i. The corresponding 
stability properties and behavior of kinetic critical thickness are more complicated than those of the system with 
planar spacer layers studied above, due to the rough strained/spacer interfaces leading to a stronger coupling of strain 
fields. 



VI. CONCLUSIONS 

We have studied the early morphological evolution and stability/instability properties of coherent multilayer struc- 
tures with the stacking of layers of different materials. Generally, the diffusion in the bulk film is assumed to be 
negligible, thus the phenomena observed in the whole multilayer films, such as the lateral composition modulations 
in tensile/compressive multilayer and the vertically correlated arrays of islands in stained/spacer system, are in fact 
surface-induced and arise through the stacking of layers with the coupling of strain fields between surface and un- 
derlying layers. We have proposed a recursive procedure to directly calculate the clastic fields of a multilayer system 
to first order of perturbations. Then, based on the surface diffusion processes, the evolution equation of surface 
morphological perturbation is expressed through the elastic and morphological properties of all the layers underneath, 
including the misfit strains and buried interface perturbations. From the exact solution of this first-order equation, 
which is in terms of material parameters and growth conditions for all the deposited layers, we have determined the 
properties of stress-driven niorphologic;al instability for both the tensile/compressive and strained/spacer multilayers. 

In this paper, we have only applied our results of early morphological evolution, derived in sections II and III, to 
two kinds of periodic multilayer structures, but more complicated structures can be studied similarly. Note that for 
multilayers with very thin layer thickness, i.e., for short-period superlattices, although our results can qualitatively 
explain and predict some experimental observations, more realistic factors such as the wetting effect should be included 
in our model to quantitatively describe the system properties. Also, studies beyond the linear analysis are needed 
to determine the details of morphological modulations or islands evolution during the multilayer growth. All these 
theoretical results can be useful in the fabrication of high-quality and novel multilayer structures. 

Acknowledgments 

This work was supported by the NSERC of Canada. 



18 



APPENDIX: RESULTS FOR FIRST AND SECOND LAYER PERTURBATIONS 

The morphological perturbation hi for first deposited layer 1 is the same as that derived in the case of single-layer 
film growth. For the planar film-substrate interface (C = 0), it is 

/ii(g,t) = e'^^*/iO(g), (A.l) 

with the assumption of identical elastic constants. Here ai = ri(£^'ef q'^ — 7ig^) and h'^{q) is the initial perturbation. 
For the layer with fixed thickness h, Eq. (A.l) leads to 

fn{q) = e'''^^/'''h°{q). (A.2) 

The perturbation /12 for second deposited layer can be calculated from the method in Sec. II and the expression (A.2) 
for hi , with the results 

Mq, t) = \ [{qv2 + Fa (iJ'eaeig^ - 72/)) e"^^* 

qv2 + 0-2 '-^ ^ " 

+ F2£'e2(e2 - e^g^e-^"^*] e'''^'''''h\q), (A.3) 

where 0-2 = T2{E'e2q^ - 72?^), and 

h2{q)=h2{q,t = h/v2) (A.4) 

when second layer has fixed thickness l2- 

In some real growth systems, ^^'^^ the substrate surface is nonplanar, corresponding to the case of rough (( 7^ 0) 
film-substrate interface. The clastic state of both film and substrate is altered by this rough interface, with the 
detailed results shown in Ref. 28. Here we use these elastic results to obtain the expression of first layer perturbation 
for single-component film: 

hi{q, t) = [{qvi - ri7ig4) e-i* + TiE'eWe-'i--'] ((q), (A.5) 

qvi+ai 

where we have assumed the initial condition hi{q,t = 0) = C{q)- Note that the contribution of nonzero interface 
morphology ( is reflected in the second term of r.h.s. of Eq. (A.5), which is exponentially damped as the film grows. 
Therefore, the occurrence of instability is still determined by the real part of perturbation growth rate ai , as in the 
case of flat fllm-substrate interface [see Eq. (A.l)],^^"^^'^^'^* and the introducing of nonzero C does not influence the 
stability/instability condition for single-layer film growth. On the other hand, Eq. (A.5) shows that compared to the 
case of flat interface C = 0, the kinetic critical thickness for the onset of 3D instability^^ (see also Sec. V for the 
definition) decreases for rough interface (^ ^ 0). Also, the nonzero substrate roughness ( can affect the detailed film 
morphology, as found in experiments. 

The second layer surface perturbation can then be calculated for rough film-substrate interface, with the result: 

h2{q,t) = ^-—l l^\{qvi-Ti^iq^)e^^'^/--+TiE'elq''e-'''^ 

qv2 + CT2 (qvi +ai ' 

X {{qv2 + F2 {E't2eiq^ - 729^)) e"^* + F2^'e2(e2 - ei)q^e-'i^-' 



+ r2i;Ve2ei(e-'"'=*-e'^^*)e-9'^| 



C(g). (A.6) 



Similar to Eqs. (A.2) and (A.4), the surface morphological perturbations of first and second layer with fixed layer 
thicknesses, i.e., /ii(g') and h2{q), can be obtained by substituting t = li/vi and t = I2/V2 into Eqs. (A.5) and (A.6), 
respectively. 



* Electronic address: huang@csit.fsu.edu; Present address: School of Computational Science and Information Technology, 

Florida State University, Tallahassee, Florida 32306-4120 
^ Electronic address: desai@physics.utoronto.ca 



19 



^ V. A. Shchukiri and D. Bimberg, Rev. Mod. Phys. 71, 1125 (1999). 

^ A. G. Norman, S. P. Ahrcnkiel, H. Moutinho, M. M. Al-,]assim, A. Mascarcnhas, ,J. Mirccki Millunchick, S. R. Lee, R. D. 

Twcstcn, D. M. FoUstacdt, J. L. Reno, and E. D. Jones. Appl. Phys. Lett. 73, 1844 (1998). 
^ S. P. Ahrenkiel, A. G. Norman, M. M. Al-Jassim, A. Mascarenhas, J. Mirecki-Millunchick, R. D. Twesten, S. R. Lee, D. M. 

FoUstaedt, and E. D. Jones, J. AppL Phys. 84, 6088 (1998). 
" R. D. Twesten, D. M. Follstaedt, S. R. Lee, E. D. Jones, J. L. Reno, J. Mirecki Millunchick, A. G. Norman, S. P. Ahrenkiel, 

and A. Mascarenhas, Phys. Rev. B 60, 13 619 (1999). 
^ A. G. Norman, H. R. Moutinho, M. M. Al-Jassim, A. Mascarenhas, S. R. Lee, J.Mirecki Millunchick, J. L. Reno, and D. M. 

FoUstacdt, Appl. Phys. Lett, (to be published). 

A. Ponchct, A. Rocher, J.-Y. Emery, C. Starck, and L. Goldstein, J. Appl. Phys. 74, 3778 (1993). 
^ A. Ponchet, A. Rocher, A. Ougazzaden, and A. Mircea, J. Appl. Phys. 75, 7881 (1994). 

* A. Ponchet, A. LeCorre, A. Godefroy, S. Salaiin, and A. Poudoulec, J. Cryst. Growth 153, 71 (1995). 
^ J. Tersoff, C. Teichert, and M. G. Lagally, Phys. Rev. Lett. 76, 1675 (1996). 

^° C. Teichert, M. G. Lagally, L. J. Peticolas, J. C. Bean, and J. Tersoff, Phys. Rev. B 53, 16 334 (1996). 
^1 Q. Xio, A. Madhukar, P. Chen, and N. P. Kobayashi, Phys. Rev. Lett. 75, 2542 (1995). 

1^ J. C. Gonzalez, F. M. Matinaga, W. N. Rodrigues, M. V. B. Moroira, A. G. de Oliveira, M. I. N. da Silva, J. M. G. Vilela, 
M. S. Andrade, D. Ugartc, and P. C. Silva, Appl. Phys. Lett. 76, 3400 (2000). 
Y. Nakata, Y. Sugiyama, T. Futatsugi, N. Yokoyama, J. Cryst. Growth 175/176, 713 (1997). 
O. G. Schmidt, O. Kienzle, Y. Hao, K. Eberl, and F. Ernst, Appl. Phys. Lett. 74, 1272 (1999). 

V. Le Thanh, V. Yam, P. Boucaud, F. Fortuna, C. Ulysse, D. Bouchier, L. Vervoort, and J. M. Lourtioz, Phys. Rev. B 60, 

5851 (1999). 

^® V. Le Thanh, V. Yam, P. Boucaud, Y. Zheng, and D. Bouchier, Thin Solid Films 369, 43 (2000). 
" F. Liu, S. E. Davenport, H. M. Evans and M. G. Lagally, Phys. Rev. Lett. 82, 2528 (1999). 
H. Lafontaine, N. L. Rowell, and S. Janz, Appl. Phys. Lett. 72, 2430 (1998). 

R. J. Asaro and W. A. Tiller, Metall. Trans. 3, 1789 (1972); M. A. Grinfeld, Sov. Phys. Dokl. 31, 831 (1987); D. J. Srolovitz, 
Acta Metall. 37, 621 (1989). 

^° B. J. Spencer, P. W. Voorhees, and S. H. Davis, Phys. Rev. Lett. 67, 3696 (1991); J. Appl. Phys. 73, 4955 (1993). 
^1 J. E. Guyor and P. W. Voorhees, Phys. Rev. Lett. 74, 4031 (1995); Phys. Rev. B 54, 11710 (1996). 

B. J. Spencer, P. W. Voorhees, and J. Tersoff, Appl. Phys. Lett. 76, 3022 (2000); Phys. Rev. B 64, 235318 (2001). 
N. Sridhar, J. M. Rickman, and D. J. Srolovitz, J. Appl. Phys. 82, 4852 (1997). 

2" L. E. Shilkrot, D. J. Srolovitz, and J. Tersoff, Appl. Phys. Lett. 77, 304 (2000); Phys. Rev. B 62, 8397 (2000). 

J. D. Eshclby, Proc. Roy. Soc. London A 241, 376 (1957). 
^® A. A. Maradudin and R. F. Wallis, Surf. Sci. 91, 423 (1980). 

F. Leonard and R. C. Desai, Phys. Rev. B 57, 4805 (1998); Appl. Phys. Lett. 74, 40 (1999). 
2* Z. F. Huang and R. C. Desai, Phys. Rev. B 65, 205419 (2002). 

B. J. Spencer and J. Tersoff, Phys. Rev. Lett. 79, 4858 (1997); J. Tersoff, Phys. Rev. B 43, 9377 (1991). 
^° P. Sutter and M. G. Lagally, Phys. Rev. Lett. 84, 4637 (2000). 
3^ R. M. Tromp, F. M. Ross and M. C. Renter, Phys. Rev. Lett. 84, 4641 (2000). 

H. Sunamura, N. Usami, Y. Shiraki, and S. Fukatsu, Appl. Phys. Lett. 66, 3024 (1995). 

H. Lafontaine, B. F. Mason, S. J. Rolfe, D. D. Perovic, and B. Bahierathan, J. Vac. Sci. Technol. B 16, 599 (1998). 
A. Ballestad, B. J. Ruck, M. Adamcyk, T. Pinnington, and T. Tiedje, Phys. Rev. Lett. 86, 2377 (2001). 



20 




1 

Ua,o|k 
(substrate) 

(a) 



to first order 




o(=0 



•■a, kk 



Ua,k-lk 



Ua, l|k 



Ua,o|k 
(b) 



Ik 
-hk-2 



Ua,kk 



s.hk-1 



identical 
elastic constants 



-ho 



Ua,k-l|k 



z«0 



(c) 



FIG. 1: Schematic diagram of the simplification procedure for determining the elastic state of the multilayer system, (a) The 
original fe-layer structure that we study, with growing and undulating surface (located at z = hk{x,y,t)) and with nonplanar 
and frozen interfaces (located at z = hi{x,y)) between buried layers i = k — 1, fc — 2, 1, and (substrate). The elastic 
displacement field of each layer is indicated as Ua,i\k with a = x, y, z. (b) The efi^ective system of which the elastic state 
is identical to that of system (a) to first order of perturbations. Note the planar interfaces between layers i and i — 1 with 
< i < fc — 1, which arc considered to be in the effective clastic media of new-state fields u2,i\k- (c) The simplified system that 
is equivalent to system (b) when the elastic constants of all the layers and substrates are assumed to be identical. Then we 
have only one layer on a semi-infinite effective medium, with a rough surface and single nonplanar interface. In both (b) and 
(c), top free surface obeys the boundary condition (16) and the k/{k — 1) interface fulfills the special boundary condition (17). 
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FIG. 2: Stability diagram of rescaled layer thickness Ia versus deposition rate uj, with parameters e*B = —1, Ib/Ia = !> 
Vb/va = !> and 7^ = 1 for the growing multilayer system. Different relative surface mobilities T*b = 0.001 (dashed curve) and 
0.01 (solid curve) are used, as indicated in the figure. Region marked as "S" is effectively stable. Inset: — 1 (solid) and 100 
(dotted). 
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FIG. 3: Stability diagram of l\ versus v\, with parameters T*g = 0.01, Vg/v\ = 1, and — 1. (a) e*g — —1, with different 
values of Ib/Ia- 1 (solid), 2 (dashed), and 3 (dotted); (b) Fixed value of 1%/Ia = Ij but different e% = —1 (solid), —2 (dashed), 
and —3 (dotted). 
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FIG. 4: Stability diagrams of Ib/Ia versus v'X, with parameters eg = —1, Vg/v'^ = 1, and ^g = 1. (a) Relative surface 
mobilities r*g = 0.01 (solid) and 0.001 (dashed), with + = 0.9 and 0.3; (b) = 100 (solid) and 1 (dashed), with 
Ia + Ib= 0.9. Different types of asymmetry can be seen for different values of r%. 
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FIG. 5: Stability diagram of misfit e*g versus v\, witfi = = 0.5, vb/va = 1; 7b = !> a-nd different — fOO (a) and O.Of 
(b), sfiowing different typos of asymmetry. Panels (c) and (d) are the replots of (a) and (b) as —1/cb versus v^, in order to 
facilitate comparison respectively with Fig. 6 and its inset, where the strain-balanced condition implies Ib/Ia = — l/es- 




FIG. 6: Stability diagram of Is /I a versus Va in the strain-balanced condition, with Vg/v^ = 1 and 7^ = 1. For relative 
mobility FJ, = 0.01, 1*a+Ib is fixed at 0.9 (solid) and 0.3 (dashed). Inset: F^ = 100 (solid) and 1 (dotted), with + = 0.9. 
Compared to Figs. 4 and 5, here the stability property shows single type of asymmetry for different Fg. 
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FIG. 7: Stability diagrams for different relative surface tension 7I; = 0.1 (dashed), 1 (solid), and 10 (dotted), with Vg/v^ = 1- 
In (a), v'X is plotted as a function of Z^, with r% = 0.01, e*g — —1, 1*b/1*a = 1, and v*g/v\ = 1; while in (b), the system is in the 
strain-balanced condition, and we use different V*g [0.01 and 1 (inset)] in the stability diagram of 1*b/1*a versus v\, with fixed 

l\+l*g= 0.9. 
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FIG. 8: Maximum perturbation amplitude A'^'^^ versus time t, for successive strained layer numbers k = 1,3, ...,11 (from 
bottom to top). Note the semilog plot. The material parameters used are analogous to those of Ge/Si(001) growth system, 
with Ia=4: ml, Ib = 9 nm, va = 1 ML/min, and T = 550 °C. 
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FIG. 9: Kinetic critical thickness hc\k (rcscalcd by that of single layer hc\k=i) as a function of the layer mimbcr fc, with the 
Ge/Si(001) growth parameters: Ia ~ 4: ML, iia = 1 ML/min, T = 550 °C, and different spacer layer thickness Ib ~ 9.5 nm 
(square), 8.5 nm (triangle), and 8 nm (circle). Inset: hc\k/hc\k=i as a function of Ib, for A; = 3, 5, 7, and 11. 
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FIG. 10: Rescaled kinetic critical thickness hc\k/hc\k=i as a function of average strained layer thickness Ia, for successive 
strained layer numbers. Ib = 8.5 nm, va = 1 ML/min, and T = 550 °C. 
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FIG. 11: Kinetic critical thickness hc\k as a function of (a) misfit eA and (b) deposition rate va, for different layer numbers 
A; = 1, 3, 5, and 11. The parameters are Za = 4 ML, Ib = 9 nm, and T = 550 °C for both (a) and (b), as well as = 1 
ML/min for (a). Note the inset in (b) where the range of va is almost ten times larger. 



